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Abstract 

We give three examples of stochastic processes in the Gelfand-Tsetlin cone in 
which each component evolves independently apart from a blocking and pushing 
interaction. These processes give rise to couplings between certain conditioned 
Markov processes, last passage times and exclusion processes. In the first two 
examples, we deduce known identities in distribution between such processes whilst 
in the third example, the components of the process cannot escape past a wall at 
the origin and we obtain a new relation. 

Keywords: Gelfand-Tsetlin cone; conditioned Markov process; exclusion pro- 
cess; last passage percolation; random matrices 

AMS 2000 Subject Classification: Primary 60J25; Secondary: 60C05. 

Submitted to EJP on May 2, 2009, final version accepted July 16, 2009. 



1 Introduction 

In pQ , the authors Baik, Deift and Johansson show that suitably rescaled, the law of the 
longest increasing subsequence of a uniformly chosen random permutation of {1,2,..., n} 
converges, as n tends to infinity, to that of the Tracy- Widom distribution. The latter, 
first identified in [2E] , describes the typical fluctuations of the largest eigenvalue of a large 
random Hermitian matrix from the Gaussian unitary ensemble (see [11] for a definition). 
This somewhat surprising discovery has been followed by much research which has shown 
that the Tracy- Widom distribution also occurs as a limiting law in various other models 
such as last passage percolation [121 EES], exclusion processes [25], random tilings [T4"l [T3] 
and polynuclear growth [Tj)J [22] • See also the survey [IB] . 

Eigenvalues of random matrices are closely related to multi-dimensional random walks 
whose components are conditioned not to collide. In particular, both fall into a class 
of processes with determinantal correlation structure and exhibit pairwise repulsion at a 
distance. On the other hand, models such as the exclusion process are defined by local 
"hard edged" interactions rather than particles repelling each other remotely. This paper 
is concerned with showing how it is possible to connect these two types of model by 
coupling processes of one class with processes from the other. 
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In common with previous works in this area, we realise these couplings via the construc- 
tion of a stochastic process in the Gelfand-Tsetlin cone 

K n = {(a; 1 , x 2 , . . . , x n ) E R 1 x R 2 x . . . x W 1 : x* +1 < x\ < x^}. 

A configuration (x 1 , . . . , x n ) G K n is called a Gelfand-Tsetlin pattern and we may repre- 
sent the interlacing conditions x\ +1 < x\ < a^i 1 diagrammatically as follows. 



x- t 



.1, ■-) 



X- 



.?•: 



x 



n-l 



x:; 



Suitable processes in the Gelfand-Tsetlin cone appear naturally in several settings, for 
example the particle process associated with a random domino tiling of the Aztec diamond 
[T5] and the eigenvalues of a GUE matrix and its minors [2]. In other cases, the process 
in K n is not evident at first sight and must be constructed, for example see the recent 
studies of asymmetric simple exclusion processes jU [6] . 

Most frequently, dynamics for the process in K n are constructed using a combinatorial 
procedure known as the Robinson-Schensted-Knuth algorithm (see O'Connell [20] )• With 
RSK dynamics, the n(n + l)/2 components of the process are driven by a noise with only 
n degrees of freedom, leading to strong correlations between components. 

In this paper we consider some alternative dynamics in which every component of the 
process evolves independently except for certain blocking and pushing interactions that 
ensures the process stays in WL n . This approach yields a new relation between an exclusion 
type process constrained by an impenetrable wall and a multi-dimensional random walk 
with components conditioned to neither become disordered nor jump over the wall. Dy- 
namics of this type have previously been considered by Warren [29] for Brownian particles 
(see also Toth and Veto [2Z]), by Nordemstam in the context of shuffling domino tilings 
of the Aztec diamond [TH] , by Borodin and Ferrari in the context of surface growth mod- 
els [3]. Analogous dynamics have also previously been studied in the context of growth 
models where they are known as Gates and Westcott dynamics, see Prahofer and Spohn 
[22] for example. 



2 Description of dynamics and results 



From here on, we work exclusively with Gelfand-Tsetlin patterns with integer valued 
components and hence modify our definition of K n to 

K n = {{x 1 , x 2 , . . . , x n ) 6 Z 1 x Z 2 x . . . x Z" : x\ +1 < x\ < x k +l). 
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2.1 Poisson case 



Our first example consists of a continuous time K n valued Markov process (X(i);t > 0) 
that determines the positions of n(n + l)/2 interlaced particles on the integer lattice Z 
at time t. The stochastic evolution of the pattern X is as follows. 

Fix a vector of rates q £ (0, oo) ra and identify each particle with its corresponding com- 
ponent in X. The particle X\ jumps rightwards at rate q\ > 0, i.e. after an exponentially 
distributed waiting time of mean q± . The two particles, X\,X\ corresponding to the 
second row of the pattern each jump rightwards at rate q2 independently of X\ and each 
other unless either 



• %\{t) = X\(t), in which case any rightward jump of X\ is suppressed {blocked), or 

• X|(£) = X\(t), in which case X\ wm be forced to jump {pushed) if X\ jumps. 

In general, for k > 1 and 1 < j < k, each particle X^ attempts to jump rightwards at 
rate q^, and will succeed in doing so unless it is blocked by particle X^ -1 . Particle X k k 
can always jump rightwards at rate q^ without impediment. In addition, if X^ -1 = X^ +1 , 
particle X^ +1 is pushed to the right when X^ -1 jumps. This blocking and pushing ensures 
that X(t) remains in K n for every t > 0. We will show that for certain initial conditions on 
X(0), the bottom layer of the pattern, (X n (t); t > 0), is distributed as a mult i- dimensional 
random walk with components conditioned not to become disordered (Theorem 12. ip . 

To describe the result more precisely, recall that for z £ W n = {z £ Z n : Z\ < z<± < . . . < 
z n }, the Schur function S z : W 1 — > R can be defined (see for example [10J ) as a sum of 
geometrically weighted patterns, 

S z ( qi ,...,q n )= w ^ x )- (2- 1 ) 

x£K n (z) 

The sum is over K n (z) = {x £ IK„ : x n = z}, the set of all Gelfand-Tsetlin patterns 
x = (x 1 , . . . , x n ) £ K„ with bottom row x n equal to z and the geometric weight function 
is 

n 

w \x) = [[ql 

i=l 

where \z\ = Ylt=i z i f° r z ^ ^ an d = 0- 

This definition gives a natural probability mass function on patterns x £ K n (z), 

«.<*) = ||. (2.2) 



Suppose that {Z(t); t > 0) is an n-dimensional random walk in which component i 
is, independently of the other components, a Poisson counting process of rate g,. The 
function h : W n — > M defined by 

Mx) = gr Xl ---^"^(9). (2.3) 
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is harmonic for Z killed at the first instant it leaves W n (see [T7j for example). Hence, h 
may be used to define a new process, Z\ with conservative Q-matrix on W n defined by 

N h(x + ei) S x+ei (q) , , . , ^ 
Q(x,x + ei) =qi \ - = - , '\ 1 < i < n, x G W n , 

where e; is the standard basis vector, and the other off diagonal rates in Q are zero. 

This Doob /i-transform, Z', may be interpretted as a version of Z conditioned not to 
leave W n and is closely related to the Charlier orthogonal polynomial ensemble (again 
see PH). 

In section |3] we prove the following result, obtained independently by Borodin and Ferrari 
by another method in [3]. 

Theorem 2.1. If (£(£); t > 0) has initial distribution M z (-) for some z 6 W n t/ien 
(X n (t); t > 0) is distributed as an n dimensional Markov process with conservative Q- 
matrix 

Q(x, x + ei ) = lix+etew*], 1 < i < n, x G W n . 

and all other off diagonal entries set to zero, started from z. 



Note that from structure of the initial distribution and the construction of X, this theorem 
implies that in fact every row of the pattern is distributed as a conditioned Markov process 
of appropriate dimension and rates. 

Theorem 12.11 readily yields a coupling of the type discussed in the introduction - the 
(shifted) left hand edge {X\(t),X\(t) — 1, . . . ,X"(t) — n + 1; t > 0) of X has the same 
"hard edged" interactions as an asymmetric exclusion process (the particle with position 
Xi(t) — k + 1, 1 < k < n takes unit jumps rightwards at rate but is barred from 
occupying the same site as any particle to its right). However, Theorem 12.11 implies that 
(£"(t); t > 0) has the same law as (Z\(t); t > 0), the first component of the random walk 
Z conditioned to stay in W n , when started from Z^(0) = z. Further we observe that 
when z — (0, . . . , 0), M z is concentrated on the origin and a version of the left hand edge 
can be constructed from the paths of Z via £\(t) = Z\(t) and 

Xl +1 (t) = Z k+1 (t)+ inf (Xt(s) - Z k+1 (s)) , l<k<n. 



Iterating this expression and appealing to Theorem 12.1} 

(Z\(t); t > 0) t» M ^ ± (Z,(t.) - Z,(t,^) ; t > ) . (2.4) 



This identity was previously derived by O'Connell and Yor in [21] using a construction 
based on the Robinson-Schensted-Knuth correspondence. 
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2.2 Geometric jumps 



For our second example we consider a discrete time process (£(£); t G Z + ) (where 
Z + is the set of non-negative integers) in K n in which components make independent 
geometrically distributed jumps perturbed by interactions that maintain the interlacing 
constraints. 

Let q be a fixed vector in (0, l) n and update the pattern at time t beginning with the top 
particle by setting X\(t + 1) = X\(t) + £, where £ is a geometric random variable with 
mean (1 — q\)/q\. That is, the top most particle always takes geometrically distributed 
jumps rightwards without experiencing pushing or blocking. 

Suppose rows 1 through k — 1 have been updated for some k > 1 and we wish to update 
the position of the particle corresponding to the j th component of the k th row in the 
pattern, If 3L*z\(t + 1) > %j{t), then 3£*(t) is pushed to an intermediate position 

£j(t) = 3eJZi(* + !)> while if ^J-i(* + !) < X j (*)> no Pushing occurs and £j(t) = £)(t). 




£?(*„) 3£?(* + 1) 3^(<o) *i(t ) ^(*o + l) 



Figure 1: Example of blocking and pushing 

The particle Xh then attempts to make a rightward jump of size that is geometrically 
distributed with mean (1 — quj/qu from its intermediate position Xj(t) (so the particle is 
pushed before it attempts to jump). It always succeeds if j = k (i.e. it is the right most 
particle) while if j < k, it cannot jump past X k ~ l (t), the position of particle to the right 
of it on the row above before the update. The leftmost particle k > 1 is not subject 
to pushing by any particle, but is still blocked by the "ghost" of the particle 3t x — 1 . 

To state the result, let us write x -< x' when the inequality X\ < x[ < x 2 < . . . < x' n _ x < 
%n < x' n holds for x, x' G W 1 and suppose M z is as defined in ( 12. 2p . Then, 

Theorem 2.2. If 36(0) has initial distribution M z (-) for some z G W n then (X n (t);t G 
Z + ) is distributed as an n dimensional Markov process in W n with transition kernel 

p{x,x') = f[(l - %)§44 t[^x>], x,x' G W n , 

i=l x ^> 

beginning at z. 

The Markov process with transition kernel p can be described by a Doob /i-transform - 
suppose Z is now a discrete time random walk beginning at z G W n in which the k th 
component makes a geometric (qk) rightward jump at each time step, independently of 
the other components. Then the function h defined in (12.31) is harmonic for Z killed 
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at the instant that the interlacing condition Z(t) -< Z(t + 1) fails to hold (see |20j). 
The corresponding /i-transform is the discrete analogue of a process that arises from 
eigenvalues of Wishart matrices [8] . 

The right hand edge of the pattern, (3C\(t), 3L\{t), . . . ,X^(t);t G Z + ) has a simple con- 
nection to the last passage percolation model with geometric weights that may be for- 
mulated as follows. Suppose that rjk(t) are independent geometric (g^) random variables 
attached to sites in the lattice 1 < k < n, t > 1. An increasing path 7r from (1,1) to 
(t, k) is a collection of sites {(ti, ki), . . . , (tjv, &tv)}; N = t + k — 2, such that the step 
(t m +i, k m +i) ~ (tm, k m ) G {(1, 0), (0, 1)}, and we denote the set of such paths by U(t, k). 
The quantity of interest is the k- dimensional process of last passage times 

Gk(t) = max \ Vjii)) t e 

It is not difficult to confirm that [G\(t), . . . , G n (t); t G Z + ) has the same law as the 
right hand edge (3t\(t), X\{t), . . . , t G Z + ) when X has initial distribution M z , 

z = (0, . . . , 0) (i.e. 3tj.(0) = 0, 1 < k < j < n). But, a version of the right hand edge may 
be constructed from paths of Z begun at the origin so that Theorem 12.21 gives 



(Zt(t);t>0) d =M max V (Z,(z)-Z^-l));t>0 . (2.5) 

V en(i - n) (fe / 

As a consequence, Theorem 12.21 provides a new proof that such last passage percolation 
times have the same distribution as the rightmost particle in the conditioned process 
Z^ (the distribution of which, at a fixed time, is given by the Meixner ensemble - see 
Johansson [13] or [T2] ). This is a key step in obtaining the Tracy- Widom distribution in 
this setting. 

Note that the dynamics discussed above are different to those exhibited in [3] for geometric 
jumps. In particular, the particles in the process we described above are blocked by the 
position of the particle immediately above and to the right of them at the previous time 
step. 



2.3 With wall at the origin 

The final example of the paper uses the ideas introduced above to construct a continuous 
time process (£(£); t > 0) on a symplectic Gelfand-Tsetlin cone. The latter are so termed 
because they are in direct correspondence with the symplectic tableau arising from the 
representations of the symplectic group [26] . 

The space K° of integer valued symplectic Gelfand-Tsetlin patterns may be defined (see 
for example [7] or [23]) as the set of point configurations (x 1 , x 2 , . . . , x n ~ l , x n ) such that 

• x 2i ~\x 2i G Z ? + for 1 < i < LfJ and x n G zf +1)/2 if n is odd, 
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• x 2i ~ x -< x 2t for 1 < i < LfJ , 

• x 2i ■< x 2i+1 for 1 < i < L^J . 

So the all the points in a symplectic pattern lie to the right of an impenetrable wall at 
the origin, represented diagrammatically below. 



.2.' -I 



.if. I 



X5 
2 



, .4 

j; 2 



.1" 



In the vein of previous sections, we construct a process (X(t),t > 0) in K° in which 
only one particle jumps of its own volition at each instant and a blocking and pushing 
interaction maintains the interlacing constraints. 

Fix q G (0, l) n . The top particle X\ jumps right at rate q\ and left at rate , apart 
from at origin where its left jumps are suppressed. The second row also only has one 
particle, X\, which jumps rightwards at rate gf 1 and leftwards at rate q± (notice rates 
are reversed), except at instances when X\(t) = X\{t). In the latter case, it is pushed 
rightwards if X\ jumps to the right and any leftward jumps are suppressed. 

The remaining particles evolve in a similar fashion - on row 2k — 1, particles take steps 
to the right at rate qu and left at rate q^ 1 when they are not subject to the blocking 
or pushing required to keep the process in the state space, in particular X 2k ~ 1 has any 
leftward jump from the origin suppressed. On row 2k, the rates are reversed but the same 
blocking and pushing mantra applies. 

We will deduce that for appropriate initial conditions, the marginal distribution of each 
row (X k (t);t > 0) is a Markov process. The Q-matrices for the marginal processes can be 
written in terms of symplectic Schur functions, the definition of which is similar to that 
of the classic Schur function fl2.ll) - they are sums over geometrically weighted symplectic 
Gelfand-Tsetlin patterns. 

Fix k > and suppose that either n = 2k — 1 or n = 2k. Now let Wq = {z G Z fc : 
< zi < z 2 ■ ■ ■ < z k } and define K°(,z) to be the set of symplectic patterns x in K° with 
bottom row x n equal to z G Wq . The geometric weight on K° is 

and 

k 

«4fc ( x ) = [[ ( li 

i=l 



fe-1 

n 

i=l 



\x 2i - 1 \-\x 2i - 2 \ + \x 2i - 1 \-\x 2 
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using the convention that \x°\ = and empty products are equal to 1 (so w[ 9l \x) = q± '). 
Then, the symplectic Schur function Sp™ : M. k — > R, z G Wq , k > 1 is defined (see [H]) by 

Sp n z (q h ...,q k )= J2 <(*)■ ( 2 - 6 ) 

xeK°{z) 

For even n, S'p 71 gives the characters of irreducible representations of the symplectic 
group Spin) [26]. For odd n, Sp n was introduced by Proctor [23] and can interpretted as 
the character of the irreducible representations of a group that interpolates between the 
classical groups Spin) and Sp(n + 1) [18J. 

Define the Q-matrix Q n : Wq x Wq — > M. as follows. For x G Wq and x ± G Wq , some 
1 < i < fc, 

Q n (x,x±e,) = ^M. (2.7) 
All other off diagonal entries vanish and the diagonals are given by 

k-l 

- Qtk-ifa x) = ^2 (ft + QT 1 ) + 4k 1 l[«ci>o] ( 2 -8) 

8=1 

and 

-Q 2fe M = ^(ft + ft r1 )- ( 2 - 9 ) 

i=l 

A corollary of the intertwinings we prove in sections 15. 1 1 and l5\2l is that Q n is conservative. 
For z G W fc define M™ : K° n -> [0, 1] by 

<0) 



M z "(x) 



Then the definitions of the symplectic Schur functions imply that M^(-) gives a proba- 
bility distribution on patterns in K^{z). 

From these ingredients we obtain 

Theorem 2.3. Suppose X has initial distribution given by M"(-) ; then {X n (t);t > 0) is 
distributed as a Markov process with Q -matrix Q n , started from z. 

The relevance of this theorem to the discussion in the introduction may again be seen 
by examining the evolution of the right hand edge of X. Suppose we have a system of n 
particles with positions (3E\(t), j£f (t) + 1, 3t|(*) + 2, . . . , 3£?/ n+1 w 2 , (t) + n — 1; t > 0). 

Particle i > 1 attempts to jump rightwards at rate 7« = ?(i+i)/2 if i is odd or 7^ = qZ 2 
if i is even and leftwards at rate 7^ . An attempted left jump succeeds only if the 
destination site is vacant, otherwise it is suppressed. A rightward jump always succeeds, 
and, any particle occupying the destination site is pushed rightwards. A particle being 
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pushed rightwards also pushes any particle standing in its way, so a rightward jump by 
a particle could cause many particles to be pushed. So far we have essentially described 
the dynamics of the "PushASEP" process introduced in [I]. Our process differs by the 
presence of a wall: the leftmost particle (identified with X\) is modified so that any 
leftward jump at the origin suppressed. Also, the particle rates are restricted in that for 
odd i, the jump rates of particle % and i + 1 are inverses of each other (which is not the 
case in jl]). 

As in the previous examples, the bottom row (X n (t);t > 0) may be realised as a Doob 
/i-transform and we deduce identities analogous to (12.41) and (12.51) . For simplicity, we 
shall only consider the case that n = 2k. The case of odd n can be treated with similar 
arguments but it is complicated slightly due to the non-standard behaviour of X™ at the 
wall. 

Let Z be a fc- dimensional random walk in which the i th component jumps rightwards at 
rate q^ 1 and leftwards at rate g^. It is readily seen that Q 2 k is the Q-matrix of Z\ the 
/t-transform of Z killed on leaving Wq under harmonic functions 

h 2k (x) = q^q^ . . . q; Xk Sp 2 x k (q), x G W fc . 

Theorem 12.31 shows that yZl(t);t > Oj has the same law as (Xl k (t);t > 0) when X is 
initially distributed according to M^ k and Z(0) = z 6 Wq. 

But if z — (0, . . . , 0), a process with the same law as the right hand edge of X can be con- 
structed from the paths of Z and a random walk Z that is independent of, but identically 
distributed to —Z. The resulting identity in distribution can be stated succinctly in terms 
of the 2 k- dimensional random walk Z(t) = (Zi(t), Zi(t), Z 2 (t), Z 2 (t), . . . , Z k (t), Z k (t)) as 
follows 

(2k 
sup V(Z,(t m )-Z i (t i ));t>0 
o<h<...<t 2k+1 =t ~ 

The Brownian analogue of this result will be considered in [5]. 



3 Proof of Theorem 2.1 



Let (X(t); t > 0) be the process on K n satisfying the dynamics described in section [2TT1 It 
is clear from this description that the law of {X n (t); t > 0} is conditionally independent 
of t > 0, j < n — 1} given {3£ n-1 (t); t > 0}. That is, the dynamics of the particle 

in row n depend on the evolution of particles in the rows above only through the particles 
in row n — 1. Hence the theorem may be proven inductively by studying only the bottom 
and penultimate layers of the pattern. 

To this end, we assume for induction that the conclusion of 12. II holds. Then, when X(0) is 
distributed according to M z (-), the bottom layer (£ n (t); t > 0) is Markovian and evolves 
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according to the conservative Q-matrix Qx defined via 

Qx(x,x + e*) = a+e ' . l^+eisw*], 1 < i < n, x G W\ 
and all other off diagonal entries set to zero. 

We will define a Markov process (X(t), Y(t); t > 0) on W n ' n+1 = {(x,y) G W n x W n+1 : 
x d± y} (recall x t< y means that yi < x\ < y%+i, 1 < % < n) in which X evolves 
according to while K evolves independently of X apart from the blocking and pushing 
interaction. One should think of (X, Y) as the penultimate and bottom layer of our 
construction in K n+1 . So, to complete the induction step it is sufficient to show that 
marginally Y is Markovian and evolves according to 

Qy\V, y + ei) = Qi ft- , 

Kv) 

for y G W n+1 , where h is given by 

%) = C* 1 • • • q~li +1 S y (qu • • • , 9n+l), 2/ G 

for some g n+ i > 0, and all other off diagonal entries vanish. The diagonal entries are 
given by 

n+1 



Qy(y,y) = -J^ft- 



»=i 



Appropriate dynamics for (X, y) are specified by the conservative Q-matrix A with off 
diagonal entries given by 



A((x,y),(x',y')) 



Qx(x,x + e { ), (x',y') = (x + e h y), x { < y i+1 , 

Qx(x, x + e^, (x', y') = {x + e h y + e i+1 ), x { = y i+1 , 

Qn+i, (x',y') = (x,y + ej), 

otherwise 



for (x, y), (x', y') G VV n, " +1 and l<i<n, l<j<n + l. The diagonal entry 
-A((x',y'), (x',y')) is given by 



n+1 



Y ft + ?n+l lW<<] +9n+l = Yl ft + ft +1 • ( 3 - X ) 

i=l i=l i=l 8=1 

Now, as an immediate consequence of the definition of the Schur function in (12. ip . we 
have 

S z (q 1 ,q 2 ,...,q n ,qn+i) = Y w9 ( x ) 

xeK„(z) 

n 

y <inti A y 11^' 

z'£W n :z'<z xESL n (z') 3=1 

z'<z 



\x.i\-\x?- x \ 



10 



So the marginal distribution of the penultimate row of particles under the initial distribu- 
tion defined in (12.21) is given by m(-, y) where y G W n+1 is fixed and m : W n,n+1 — > [0, 1] 
is defined by 

/ \ \y\-\x\Sx\Q) 

m(x,y) = qn+1 

where q = (q lt q 2 , . . . , g n , g n+ i). 
Furthermore, 

Hv, W, y')) = m(x', y') l[ y >=y\ . (3.2) 

defines a Markov kernel from W n+1 to yy n ' n+1 . That is, for each y G W n+1 , A(y } •) defines 
a probability distribution on VV n ' n+1 . 

The heart of our proof is showing that the conservative Qy is intertwined with A via A, 

Q Y A = AA (3.3) 

From here, lemma IA.1I shows that A intertwines the corresponding transition kernels. 
That is, if (p t ; t > 0) are the transition kernels corresponding to Qy and (q t ; t > 0) those 
to A, then for y G W n+1 , G W n ' n+l and t > 0, 

p t (y,y')m(x',y') = ^2m(x,y)q t ((x,y), (x',y')), 

x<y 

An immediate consequence of this relationship is that for bounded / : W n+1 — > M, 

(a;',?/') x<y 

= Yl pt(y^y') m ( x '^y')f(y') 

(x',y>) 

= ^2pt(y,y')f(y') Y m ( x ',y') 

y' x'^y' 

= 5>(v,vO/(iO, 

s/' 

where E y is the expectation operator corresponding to the measure under which (X, Y) 
has initial distribution A(y, ■). 

When < t 1 < . . . < t N and /i,...,/jv : W n+1 — > R are bounded, the preceeding 
argument generalises and the intertwining shows that 

EVi(Y(h))...f N (Y(t N ))} = Pt{y,y l )vt{y\y 2 )...Pt{y N ~\y N )h{y 1 )...f N {y N ).' 



This is essentially the argument of Rogers and Pitman [21] and establishes 

Theorem 3.1. Suppose (X(t),Y(t);t > 0) is a Markov process with Q-matrix A and 
initial distribution A(y, ■), for some y G W n+1 . Then Qy and A are interwined via A 
and as a consequence, (Y(t);t > 0) is distributed as a Markov process with Q-matrix Qy, 
started from y. 
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The intertwining (13.31) is equivalent to 

Q Y (y,y') =Y^^iylA((x,y),(x',y% yEW n+1 , (x', y') E W n ' n+ \ (3.4) 
/ — ' mix' y') 

x<y y ' y ' 

where the summation is over the points x in W n that interlace with y. As the particles 
can only make unit jumps rightwards, both sides of the expression vanish unless either 
y' = y or y' = y + ej, for some 1 < j < n + 1. 

We first consider the case when y = y', corresponding to the diagonal entries of Qy. The 
right hand side of the expression is 

^m(^ 
z — ' mix' y ) 

Using the definition of m, this becomes 

Y,^ yA ^A{{x,y%{x\y')). (3.5) 

x<y' 

Now, A((x, y f ), (x',y f )) is non zero for x < y' only if x = x' or x = x' — e, for some 
1 < i < n. When x = x f , —A((x,y'), (x',y')) is the rate of leaving at (x',y'), given in 
(13. ip . On the other hand if x — x' — e^, A((x,y'), (x',y r )) is the rate at which the i th 
X particle jumps rightwards (without pushing a Y particle). But, such values of x are 
included in the summation only if x — x' — -< y — y', i.e. x\ > y\. 

Combining this with (I3.5P and (13.11) and the fact that ^ x 6l ' = q n +i, we see that if 
y = y' the right hand side of (13. 4p is 

n O f \ n+1 n 

Qn+1 s Qx{X -ei,X) l[x(>y^} ~ <?< ~ 1 W<*'J " 

i=l x * ' i=l i=l 

The first summand above is 

9n+i g 7^ Qx{x -e h x) 1[^>^] = Qn+l My'i<x'i\i 

so the first and last summations above disappear and we are left with — Y^=i which 
is exactly Q Y (y',y')- 

If y 7^ y', the only other possibility is that y' = y + for some 1 < i < n + 1. Let us 
first deal with the simplest case, where i — 1, that is, y' = y + e±. The only value of x for 
which A((x, y' — ei), (x', y')) is non zero is x = x' as the first Y particle is never pushed 
by an X particle. Furthermore, y[ — 1 < y[ < x[ and so the jump of Yi is certainly not 
blocked. Hence, 

E m \ x ->y) \ / i i\\ m(x',y' — ei) .. , , , , 

- A {{x,y , (x,y)) = . > A{{x ,y -ei),(x,y)) 
x . x ^ y m ( x >y) m{x',y>) 



Qy(y' - ei,y')- 



^T^n+l 
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So in this case, (13.41) is satisfied. 



For i > 1, consider the dichotomy x' i _ l < y[ or = y[. Suppose we are in the former 
case, i.e. y' = y + e, and < It is not possible that the movement in the i th 
component of Y could have been instigated due to pushing by the (i — l) th X particle (a 
push could only have occurred if a;^_ 1 — 1 = y[ — 1). Thus, as in the i — 1 case above, 
A((x, y' — e»), (x', y')) is non zero only for x = x' and almost identical calculations verify 

AH- 

The second i > 1 subcase is that x^„ x = and y = y' — Here the only possibility is 
that the Y particle "did not jump but was pushed" , which one may confirm by noting 
that x' does not interlace with y' — when x\_ x = y[. So, the right hand side of f)3.4p is 
given by 

m{x — ei-i,y — ej) t , . . r 

— — A{{x -ei~i,y -ei),{x,y)). 

m{x, y ) 

Using the definitions of m and A, this becomes 

Q l n+T Hx '~ ei - 11 S^.M S y/ (q) S x ,(q) 



a quantity which is easily seen to equal Qy(y' — y')- 

This concludes the proof that Qy and A are intertwined via A. 



4 Proof of Theorem 2.2 



It is again sufficient to consider any pair of consecutive rows (X, Y) and construct the 
process iteratively. 

Let (X(t);t G Z + ) be an n dimensional Markov chain in W n with one step transition 
kernel 

/ /\ /x S x i (q) 

PX{X,X) = a W) ^T-jT-j- l[x^x'] ■ 

where q G (0, l) n , a(q) = YYi=i(^ ~ ft) anc ^ f° r x i x ' e K", x -< x' indicates that the 
inequality x\ < x[ < Xi < . . . < x' n _ x < x n < x' n holds. 

Let (t G Z + , 1 < i < n+1) be geometric (q n +i) random variables that are independent 
of each other and of X, 

P(6(*0 = j) = (1 - q n +i)qi +1 , J = 0, 1, 2, ... . 

Define a process (Y(t)]t G Z + ) in W n+1 in terms of X using the recursion 

Yi(i + 1) = min(Yi(t) + + l),X a (t)), 

K n+1 (t + 1) = max(y n+1 (t), X n (t + 1)) + + 1), 

Yj{t + 1) = min(max(y j (t),X j _i(t + 1)) + + l),Xj(t)), 
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for 2 < j < n. 

The recursion encodes the blocking and pushing mechanism, maintaining the initial in- 
terlacing relationship, so X(t) -< Y(t) for each t. 

We will prove that if A is as defined in (13.21) then 

Theorem 4.1. If (X,Y) is initially distributed according to A(y, ■), y G W n+1 ; and then 
evolves according to the recursion above, the marginal process (Y(t);t > 0) is distributed 
as an n + 1 dimensional Markov process with transition kernel 

/a / \ S y '(qi, . . . , q n +i) ^ 

by{Ql, • • • > qn+l) 

started from y . 

Our strategy, again, is to prove that A interwines the corresponding transition prob- 
abilities. Suppose (x,y), (x',y r ) G W n,?1+1 , x -< x' and y -< y'. Let us write down 
q((x,y), (x',y')), the one step transition probabilities for (X, Y). Firstly note that 

q((x, y), (x, y')) = r(y, x' , x, y)p x (x, x 1 ). 

where 

r{y\ x', x, y) = P(F(1) = y'\X{l) = x', X(0) = x, Y(0) = y). 

Using the definition of Y, r can be conveniently expressed in terms of the "blocking" and 
"pushing" factors b, c : Z 2 — > M + 



b(u,v) = (1 - q n+ i) t [v<u] + t[ u=v ], 

c( u , v) = q-^ l[ u>v ] . 

Then r(y', x', x, y) is equal to 

qH^Kxi, y' x ) Y[qn+iK x i>y'i) c ( x i-i,yi) fc^ 1 - Qn+i)c{x' n ,y n+l ). 



,i=2 



To prove the theorem we will need the following "integrating out" lemma. 

Lemma 4.2. Suppose v 2 ,v[,u' G Z satisfy v[ < v 2 , v[ < v! and q n+ i ^ 1. Then we have 

V'2Au' 

QnliKn, v[)c(u', v 2 ) = q~tr 2 (4-1) 



The lemma may be understood more readily by imagining that we are considering the 
n = 1 case, so that there is one "X" particle nestled between two U Y" particles. We may 
fix the initial and final positions of the U Y" particles (v and v' in the lemma above) and 
also the final position of the "X" particle (u in the lemma) - it is the starting location 
of the X particle that we are integrating out. The summation is over the possible values 
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that the X particle may have started from. It must be at least equal to the final position 
of the left most Y particle v[, as this particle cannot overtake the X particle (see recursion 
equations above). Also, it cannot exceed either the initial position of the second Y particle 
i>2 (due to the interlacing constraint) or the final position of the X particle u' (as the 
particles may only jump rightwards). 

Proof. After using the definitions of b and c, the sum becomes 

1)2 Alt' 

^2 (( X ~~ qn+1 "> 1 K<«1 + (<?n+l l[«'<t*] +Qn+1 l[u'>wa]) 

Now expand the brackets in the summand and sum the terms individually. We find 



i>2 An' 




u=v'. 



V2Au' 

(1 - q n +l) Cff " 1[«'>d 2 ] = C+l (<7n+l - Qn+lj Mu'>v 2 ], 

u=v' 1 

v 2 Au' 

-a n -v2-u -ii _ — "2— uj t 

9n+l i[ii'<t) 2 ] — Qn+1 l[u'<v 2 ] 

u=v[ 

and 

i>2Au' 

i[w=^] 9 n+ l 1[«'>d 2 ] - <?n+l i[u'>« 2 ] • 

Summing the above expressions gives the result. □ 

The interesting thing about this scheme, as we will see in a moment, is that we may apply 
it successively from left to right when there are n particles so that the leftmost particles 
get heavier and heavier until we have reduced the problem to the n = 1 case. 

When the initial distribution is A(y, •), the joint distribution after one time step is given 
by 

n( x ', y') = m ( x ' fMfo f )> ( x '> v'))- 

x-<y 

Expanding the sum and incorporating the conditions y[ < X{ and x -< x' into the sum- 
mation indices yields 

y n +iAx' n y 2 /\x'- L 

K(x',y')= Y ■ ■ ■ Y m ( x >y)(i(( x >y)>( x ',y'))- ( 4 - 2 ) 

x n=y' n X!=y[ 
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The summand in (14.21) equals 



A=2 



xc{x n ,y n+1 )a(q)- 



s y (q) ' 

for x ^ y, x' z< y', x -< x', y -< y' and vanishes elsewhere. 

Now, one notices that we may use lemma 14.21 to iteratively evaluate the summation over 
Xi,X2, ...,x n (in that order). More concretely, first apply the lemma with v! = x[,v = 

(yi,V2),v' = (y'^y'2) to reveal that the sum Y%^ m ( x iV)Q((. x iV)i { x \v')) is e Q ual to 
y' 1 +y l 2+y3+-+yn+i+y' n+1 -x' 1 -x2-...-x n , 

Qn+i H x 2, y 2 ) 



X 



n <Fn+lK X h ViH^i-l, Vi) j c 0n) yn+ ^ a ^^fJ~^- 



v i=3 



This expression is again in a suitable form to apply lemma I4.2[ but this time with u' = 
x' 2 , v = (7/2, 2/3), v' = (7/2 5 2/3) an d summing over x 2 - Continuing in this fashion shows that 
P~2j) is equal to 

Hence we have verified the intertwining 

m(x', y')p Y (y, y) = ^2 m ( X) v)> v'))> 

and Theorem 14.11 follows from the argument of [21] discussed in the previous section. 



5 Proof of Theorem 2.3 



As in the previous two examples, we give a row by row construction. This time the 
asymmetry between odd rows and even rows means we have to specify how to iterate 
from even rows to odd rows and odd rows to even rows separately (presented below in 
15.11 and 15.21 respectively) . 

En route to proving Theorem 12.31 we need to conclude that Q n is a conservative Q-matrix 
for each n. 

This will be achieved by an inductive argument. Let H(n) denote the hypothesis that Q n 
is a conservative Q-matrix. It is easy to establish H(l), that Q\ is conservative - recall 
that for X\ > 0, Sp}s = q* 1 so 



Qi(x,x + ei) + Qi(x,x - ei) + Qi(x,x) = + l[ Xl>0 ] -qi - gf 1 l[ Xl >o], 
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a quantity equal to zero, and the off diagonal entries are clearly positive. 

Under the assumption that H(2n — 1) holds we will define a conservative Q-matrix Ao 
on Wq'" = {(x,y) G Wq x Wq : x -< y} in terms of Q^n-l anc l prove the intertwining 
relationship 

Q 2n A = AA . 

where A is a Markov kernel. Expanding the intertwining and summing both sides shows 
that ^2 X , Q2n(x,x') = 0, so we conclude that H(2n) holds as well. The step from H(2n) 
to H(2n + 1) follows a similar argument. 



5.1 Part I: Iterating from an odd row to an even row 



Suppose H(2n — 1) holds and identify Q x = Qin-\- Introduce a Q-matrix Ao on Wq : 
with off diagonal entries defined by 

Qx(x,x±ej), (x',y') = (x ± e s ,y) 
Qx(x, x - e i+1 ), (x r , y') = (x - e i+1 , y - e*), x i+1 = y t 

Qx(x, x + ej), (x', y') = (x + e j ,y + e 5 ), Xj = y$ , 

Qn\ (x',y') = (x,y±e j ) 

otherwise 



Ao{{x,y), {x',y')) 



for (x,y), (x',y') G Wq'", 1 < i < n, 1 < j < n. The diagonal entry — Ao((x, y), (x,y)) is 
given by 



n-l 



n-l 



+ Qn + q n 1 l[ Xl >0] + >_v l[j«<!Sj + i] +<?n l[j«>^]) + ?n lfe/n>Sn] +<?n 1 > (5- 1 ) 

1=1 i=l 

so under the assumption that Qx is conservative, Aq is also conservative. 
Define m : W n '" -> [0, 1] by 

\Spfr\q) 



m(x, y) = lyl 



Sp 2 v n (q) ' 



Note that the geometric factor is now q„ ' y ' instead of the usual gjf . By definition 



^ n (?) = E <(*) 

xeK° n (z) 

n—l 

E l! 2 "- 1 !-]! 2 !^"- 1 !-]^ 2 ! TT ,, 2 l* a< ~ 1 l-|3 a< l-|s a<-a 
"n "n X X * 

n-l 

= e ^ ,|Hz| e g r~ iH ^ 2| n^~ 1H * : 

= E gf'-'^ii-H?)- 
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Hence, m gives a Markov kernel A from Wq to Wq defined by 

A(y, (x',y)) = m(x',y) l [y/=y] . 

We then have 

Theorem 5.1. Assume Qm-\ is a conservative Q-matrix and (X(t),Y(t);t > 0) is a 
Markov process with Q-matrix Ao and initial distribution A(y, •) for some y G Wq . Then 
Q2n is a conservative Q-matrix and (Y(t);t > 0) is distributed as a Markov process with 
Q-matrix Q211, started from y. 

Suppose Qy = Q2n, then as usual we prove an intertwining relationship 

Q y A = AAo. 

This is equivalent to 

Q Y (y, y ') = Y^^^A ((x,y),(x',y')), yeWl {x\y') eW n ' n . (5.2) 
L — ' mix' y ) 

x<y v ' a ' 

where the sum is over x G Wq such that (x, y) G Wq 

Particles may take unit steps in either direction so we need to check the equality (j5.2p 
holds for y = y', y = y' + ej and y = y' — ej for some 1 < j < n. 

Let us first consider the case y = y'. When x = x', —Ao((x,y), (x',y')) is the rate of 
leaving (x', y') and is given by (15. ip . The only other possible values of x in the summation 
for which the summand is non-zero are x = x' ± e^, 1 < i < n. For such values (i.e. if 
(x' ±Ci,y) G Wq'"), the summand is 

m(x'±e t ,y>) gjf^^gpg^g) Spf{g) Sp^\g) 

m(x>,y>) q WHV\ S p%{q) Sp^\q) Sp^qY 

which is a rather fancy way of writing q^ 1 . But, for (x',y') G Wq'™, 

• (x' + e^y') G W n ' n onlyif^<^ 

• (x' - y') G Wq i > 1 only if x\ > y[_ x 

• (x' - e x ,y') G Wq only if ^ > 0. 



So 



Em(x,y) ,, . , . ... v-^ / _i \ 

-— — „4 ((a;, y),{x,y)) = ^ (?n +<?„ V^wl J + 



n-l 



, , m(x', y'Y 

x<y,x^x' v ' i=l 



q n t{x' n <y +Qn 1 [a;i>0] 
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On subtracting the rate of leaving — Ao((x', y'), (x', y 1 )) defined in (15.11) we find that the 
indicator functions all cancel and the right hand side of (15.21) is 



i=l 

which is equal to the left hand side. 

Next we consider the case that y' = y — ti G Wq . If i — n, the only possibility is that the 
Y particle jumped by itself. When i < n, the only possibilities are that the i th component 
of Y was pushed by the (i + l) th component of X (i.e. x = x' + ej+i) or it jumped by 
its own volition (i.e. x = x'). The former only occurs if y[ = x' i+1 , while the latter can 
only occur if y[ < x' i+1 , inducing a natural partition on the values we have to check the 
intertwining on. When y' — y — 6j, y\ < x' i+1 , i < n, or i = n, the right hand side of (15. 2\i 
is 

m(x '' y ' + <*A ffc' v' + e) W „')) - q " Ht/+ " 1 S "' rl(q) SV « (q) , - SP »" (,) 



When y' = y — e^, x' i+1 = y' { , i < n, the sum on the right hand side of the intertwining 
involves a single term, 

m(x' + e i+ i,y' + e { ) , , 

r-r-T\ Qx( x +e i+ i,x). 

m[x , y') 

Using the definitions of m and Qx shows this summand is 

ql x ' +ei+lHy ' +ei] Sp%+e i+ M) Spf{q) Sp$-\q) 



ql X ' Hy ' 1 Sp^ +ei (q) S P lr\q) S P %l i+i {qY 
Both of these quantities are equal to Qy(v' + £i, y)- 

Finally we consider the case y' — y+ei, 1 < i < n. As in the previous case, the dichotomy 
x\ = y[ and x\ < y[ divides the possible values of x in the summation into two cases, each 
of which having only one term contributing to the sum. When x\ = y' { , the i th Y particle 
must have been pushed, and 

m(x'-e u y'-e>) gjf^^gg^g) Sfflq) Sp 2 jr\q) 

m(x>, y>) W «> X > qW-W'l Spl?_ ei (q) Sp^q) Sp^q) ' 

Simplifying the expression on the right hand side by cancelling common factors in the 
numerator and denominator reveal it to be simply Qy{y' — £i,y)- 

On the other hand, when x\ < y[ the i th Y particle cannot have been pushed so the right 
hand side of the intertwining (15.21) is 

m{x '- y ' ~ ei) M( X ' v' - e) W „')) - jr'H"'-"' S P l?-\ q ) Sfffa) spin q) 

S10\\A ) y i y ) ) — \ T '\-Ui'\ C„2n /„\ a 2n-l/ 
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The proof of the intertwining relationship is concluded by noting that this is Qy{v' ' — y) 
as required. 

Now, summing both sides of the intertwining 

£ QY(y,y') m ( x> ,y') = £ ^2 m ( x ^y) A o((x,y),(x',y')) 

(x',y') {x',y') x^y 

over all pairs in (x',y') in Wq'" shows that Qy is conservative as J2 X , m(x , ,y') = 1 and 
E(,'y ) - 4 o((x,l/),(x', 2 /')) = 0. 

We then apply lemma IA.1I to recover the rest of the theorem. 



5.2 Part II: Iterating from an even row to an odd 

Suppose Qx = <?2n, Qy = <?2n+i and Ao is a conservative Q-matrix Ao on Wg' n+1 
{(x, y) E Wg x Wq +1 : x ^ y} with off diagonal entries given by 

Q x {x,x±ei), {x',y') = (x±e h y) 

Qx(x, x-ei), (x f , y') = (x - e i; y - e*), x { = y { 

Qx(x,x + ei), (x',y r ) = (x + e iy y + e i+1 ), Xi = y i+1 , 

Qn+v {x',y') = {x,y±ej) 

otherwise 



A ({x,y),(x',y')) = < 



for (x, y), (x', y') G Wq ,n+ , 1 < i < n, 1 < j < n + 1. The diagonal — Aq((x, y), (x, y)) is 
given by 

n n 

£(?« + QT 1 ) + £(<?n+l IfcKxi] +Qnll l[:y I+ i>^]) + + l[yi>o] • (5.3) 

i=l i=l 

Hence Ao is conservative if Qx is. 
From definition (12.61) we calculate 

Sp 2 z n+1 (q) = £ < +1 (g) 

^'G K 2n + l( 2 ) 

n 

E| a ,2n+l|_| a .2n| -p-j- 



| a ,2n+i|_| a ,2n| -p-j- 2 ]a; 2i - 1 1- |a: 2i |- \x 2i ~ 2 \ 



e 1 e n* 



So the function m : Wq ' n+1 — > [0, 1] given by 

,|«|-|.| Wfo) 



m(x,y) = q. 



n+1 Spl^(q) 
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induces a Markov kernel from Wq + to Wq ,n , 

A(y, (x',y')) = m(x',y') l [y/=y] . 

Our final theorem is 

Theorem 5.2. A ssume Qin i> s a conservative Q-matrix and suppose [X(t),Y(t)]t > 0) 
is a Markov process with Q-matrix Aq and initial distribution A(y, •) for some y G Wg +1 • 
Then Q2n+\ is a conservative Q-matrix and (Y(t);t > 0) distributed as a Markov 
process with Q-matrix Q2T1+1, started from y. 



The intertwining via A is equivalent to 

Q Y (y,y>) = J2 I! T^M(^y),W,y% yeWS* 1 , (x',y') e W n ' n+ \ (5.4) 

z — ' mix' y ) 

x<y y ' y ' 

where we sum over x such that (x,y) G Wq' u+1 . 

We only need to check (15. 4p holds for ?/ of the form y = y' , y = y' ± for 1 < j < n + 1 
as both sides vanish otherwise. 

Again we start with the case y' = y. When x = x', the rate of leaving —Aq((x', y'), (x', y')) 
is given by (15. 3p . The only other possible values of x for which the summand is non-zero 
are x = x' ± for 1 < % < n. For such x values satisfying (x, y) G Wq ' n+ , the definitions 
of m and Qx give 

m(x'±e h y') , ggUg) S P 2 y ? +1 (q) Sp$(q) 

m(x',y') W &l ' X) q \v'\~\*'\ Spf + \q) Spj?{g) Sp^ ei (qy 

which is equal to q^+i- But, for (x',y') G Wo'™ +1 , 

• (V + ej, y') G Wq only if x • < and 

• (x 1 - e u y') G Wq only if x- > y[. 



So, 

^ ^i^ 1 ^" 4 ^^' ^' = S fe+1 Vi<W+i] +9n+l M>Vi]) ■ 

x^y,x^x> v 5 i/ J i=1 

If we now subtract the rate of leaving (15.31) we find that at y = y' the right hand side of 
( 15. 4p is equal to 

i + <li 1 ) + q.n+i + q n li %>o] , 

which is equal to Qy(y', y')- 

The remaining cases are y = y' ± for some 1 < z < n + 1. Let us deal with y' = y — e^. 
If i = n + 1, this case corresponds to a leftward jump in the rightmost Y particle, a 
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situation that cannot arise through pushing by an X particle. If % < n + 1, then the 
jump arose by pushing if x\ — y[, while if x\ > y\ then the Y particle jumped by its own 
volition. In the case of pushing (i < n + 1, x\ — ?/•), familiar calculations show 

m(x' + e u y' + d) . , Sp 2 ? +l (q) 

MAS) QX( ' +e " X>= S^f) = QYiy - y) ' 

In the case of no pushing, i.e. % < n + 1 and x\ > y[ or i — n + 1, the summand is 

m(x',y' + ei) w , ,» Spp +1 (q) 

m(x>,y>) M{X ' V + ^ {X > y)) = SpJ^) = QAy > V } - 

Finally we consider the case y' = y + e^, 1 < % < n + 1 corresponding to a rightward 
jump in the i th Y particle. For % > 1, consider the dichotomy x' i _ 1 = y\ or x\_ x < y[, 
corresponding to the i th Y particle being pushed upwards by the (i — l) th X particle and 
a free jump respectively. The case % — 1 corresponds to the leftmost Y particle jumping 
rightwards, an event that cannot arise as a result of pushing. In the case of pushing, i.e. 
i > 1 and x\_ x = y[, the summand is equal to 

m(x — 6i—i,y — ej) ; , . . , 

— — A {{x -ei-!,y -ei),(x,y)). 

m(x', y') 

Using the definitions of Aq and m, this is 

^\-V-«-i\Si%_^{ q )Sp% + \q) „ , /N _ Spl? +1 (q) _ „ ( 

£'r M Sp^{q) SplHq) QX{X ^ } ~ Sp^q) ~ Qy ^ Vh 

as required. 

If y\ > x' i _ 1 (i > 1) or i — 1, then the i th Y particle jumped of its own accord and the 
only term in the summation is 

m(x',y'-ei) , , , , Sp 2 ? +1 (q) 

, , ,\ A {(x ,y -ei),(x,y))= - y 2n+l ,„, =Q Y {y,y)- 

This concludes the verification of the intertwining relationship and the theorem follows. 



A A lemma on intertwinings of Q-matrices 

Lemma A.l. Suppose that L and L' are uniformly bounded conservative Q-matrices on 
discrete spaces U and V that are intertwined by a Markov kernel A : U x V — > [0, 1] from 
U to V, i.e. 

LA = AL'. 

Then the transition kernels for the Markov processes with Q-matrices L and L' are also 
intertwined. 
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Note that we always use the lemma with U = W, V = W x W' where W C Z n and 
either W C Z n or W C Z n+1 . Our Markov kernel A is always such that A(u, (v, v')) > 
only if u = v', but this, of course, is not necessary for the lemma. 

Proof. The intertwining relationship LA = AV may be written 

A ( M > v)L'(v, y) = y~] L(u, u)A(u, v), u eU,v eV. 

Let (pf,t > 0) denote the transition kernels for the Markov process corresponding to Q- 
matrix L and fix uq G U. Multiplying both sides of the expanded intertwining relationship 
above by pt(u , u) and summing over u EU gives 

^2p t (u ,u) ^A(tt, v)L'(v,v) = ^2p t (u ,u) 22l(u,u)A(u,v), veV. (A.l) 

Now, let c G M be a uniform bound for the absolute values of the entries of L and V . 
Then 

\p t (u ,u)A(u,v)L'(v,v)\ < cA(u,v)p t (u ,u), 
so the double sum on the left hand side is absolutely convergent. Also, 

^ \pt(u , u)A(u, v)L(u, u)\<22 \Pt(uo, u)L(u, u) \ < 2cp t (u , u) 

u£U u£U 

and the same conclusion holds for the double sum on the right hand side. So, we may 
exchange the order of the sums on both sides to give 

^ L'{v, v) ^2p t (u ,u)A(u,v) = ^2A(u,v) ^2p t {u ,u)L(u,u), v G V. 

sev ueu ueu ueu 

Now, as (pt',t > 0) is the transition kernel corresponding to the Markov process with Q 
matrix L, it satisfies the Kolmogorov forward equation 

^_p t (u ,u) = ^2p t (u ,u)L(u,u). 

ueU 

Let us define qt{v) = J2ueu Pt( u 0j u)A(u, v) for v G V. 

We may differentiate the summation term by term in t using Fubini's theorem and the 
absolute bounds on the summands discussed above. Hence, the right hand side of (lA.ip 
is simply f t q t {v). 

Then, using the definition of qt in the left hand side of (LA. II) . we see that 

j t qt{v) = Y J L\v,v)q t (v), t > 0. (A.2) 

Now let ijp' t ] t > 0) denote the transition kernels of the Markov process with Q-matrix L', 
and 

Pt( v ) = ^(lo(v)Pt(v,v) = ^A(u ,v)p' t (v,v). 
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Then p' (v) = qo(v) for all v G V and p' t also satisfies the forward equation flA.21) in L'. 

But when the rates are uniformly bounded there is exactly one solution to the forward 
differential equation with the same boundary conditions as q t so q t (v ) = p' t (v) for all t > 
and v G V. 

By definition of qt(v ), we then have 

^2p t (u ,u)A(u,v) = ^2A(ua,v)p' t (v,v), 

and since the argument holds for arbitrary Uq G U we're done. 

□ 
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